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Abstract 

The focusing of electron flow in a symmetric p-n junction (PNJ) of graphene ribbon with different 
chiralities is studied. Considering the PNJ with the sharp interface, in a armchair ribbon, the 
electron flow emitting from (-L, 0) in n-region can always be focused perfectly at (L, 0) in p- 
region in the whole Dirac fermion regime, i.e. in whole regime Eq < t where Eq is the distance 
between Dirac-point energy and Fermi energy and t is the nearest hopping energy. For the bipolar 
ribbon with zigzag edge, however, the incoming electron flow in n-region is perfectly converged in 
p-region only in a very low energy regime with Eq < 0.05t. Moreover, for a smooth PNJ, electrons 
are backscattered near PNJ, which weakens the focusing effect. But the focusing pattern still 
remains the same as that of the sharp PNJ. In addition, quantum oscillation in charge density 
occurs due to the interference between forward and backward scattering. Finally, in the presence 
of weak perpendicular magnetic field, charge carriers are deflected in opposite directions in the 
p-region and n-region. As a result, the focusing effect is smeared. The lower energy Eq, the easier 
the focusing effect is destroyed. For the high energy Eq (e.g. Eq = 0.9i), however, the focusing 
effect can still survive in a moderate magnetic field on order of one Tesla. 

PACS numbers: 73.63.-b, 73.23.Ad, 73.40.Gk, 
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I. INTRODUCTION 



Graphene is a single layer carbon atoms packed into honeycomb lattice. From the point 
of view of its electronic properties in the low energy regime, a graphene sheet is a two- 
dimensional (2D) zero-gap semiconductor with the conical energy spectrum around Dirac 
points, the corners of the hexagonal first Brillouin zone, and its quasi-particles are formally 
described by the massless Dirac equation where the speed of light is replaced by the Fermi 
velocity of graphene.- The detailed electronic properties of graphene has been reviewed in 
Ref.j^]. Different from the usual zero-gap semiconductor in which the electrons and holes 
are normally described by separate Schrodinger equations with generally different effective 
masses, the electrons and holes in graphene are conjugately linked and described by different 
components of the same spinor wavefunction,- which means they are interconnected as 
Dirac fermions in QED. So graphene is a relativistic counterpart in the condensed-matter 
system. So far, 2D graphene has been successfully fabricated experimentally.- 1 ^ By varying 
the gate voltage^ or doping the underlying substrate,- the charge carriers of graphene can 
be easily tuned, the controllable ballistic PNJ or PNPJ are also realized experiment ally- 
Therefore intriguing phenomena exhibited in the bipolar graphene,- - — such as microwave- 
induced reflection,- specular Andreev reflection,-^ Klein tunneling,— Klein backscattering^ 
and negative refraction index effect-—, are possible to be verified experimentally. In fact, a 
direct experimental observation of Klein tunneling has been realized through a extremely 
sharp graphene PNJ.— 

It was shown that due to the Berry phase 7r, which was derived from the intersection of 
the energy bands at Dirac points,— the backscattering is absent.— This naturally leads to 
the so-called Klein tunneling^ or interband tunneling that an incident electron tunnel from 
the conduction into the valence band without backscattering. Because of the conservation 
of momentum and energy, interband tunneling through the p-n interface may resemble 
the optical refraction at the surface of meta-materials with negative refractive index.— In 
another word, the Klein paradox gives rise to the negative refraction. 19 This means that an 
interface of the symmetric PNJ perpendicular to the current flow is able to focus the electric 
current whereas a ballistic strip of p-type graphene separated by two n-type regions acts 
as a lens. These intriguing phenomena have been described in Ref . 13j , in which the Kubo 



formula was applied to the single-particle Dirac-like Hamiltonian of graphene. It means 
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that for an infinite 2D graphene system with ideal conical energy spectrum, i.e., in the very 
low energy regime, the electrons emitting from source are perfectly focused at the mirror 
symmetric point of the symmetric PNJ. For the realistic graphene system, however, one 
can not separate electrons and holes close to Dirac point due to the electron-hole puddles^ 
(about tens meV—). Of course, we can experimentally increase the density of electrons and 
holes by the gate voltage to evade from the puddle region. However, in the high density case, 
the energy spectrum deviates from the linear relation. Hence, the effect due to the non-linear 
dispersion should be examined. For this purpose, we will use the tight binding Hamiltonian 
to study the graphene based PNJ. In addition, considering the chirality of graphene ribbon, 
it is better to use the tight-binding Hamiltonian to describe the transport processes along 
different chirality directions. 

In this paper, using the tight-binding model, we carry out a theoretical study on the 
focusing effect of electron flow in the graphene ribbon with a symmetric PNJ. Due to the 
chirality of graphene, the focusing effects may be different for the zigzag ribbon and the 
armchair ribbon. Indeed, it is found that for the armchair ribbon with a sharp p-n interface 
the electron flow emitting from (-L, 0) in n-region can always be focused perfectly at (L, 0) 
in p-region for all energy E = \Ep — E p / n \ < t. For the zigzag ribbon, however, the electron 
flow is perfectly focused only in the very low regime (Eq < 0.05£). Furthermore, the perfect 
focusing in the bipolar ribbon is robust against disorders induced by the random potential. 
But the edge disorders drastically affect the perfect focusing. For a smooth PNJ, electrons 
are backscattered close to PNJ at a distance proportional to k y and interface width d— 
In this case a quantum interference between forward and backward scattering is present 
and the intensity of the focused spot is weakened, but the focusing pattern keeps almost. 
Finally, in the presence of a weak perpendicular magnetic field, the momentum k y is no 
longer conserved. Consequently particles are deflected in opposite directions in the p-region 
and n-region, which destroys the perfect focusing especially in the low energy regime. 

The rest of the paper is organized as follows. In Sec. II, the model system including 
bipolar graphene ribbon in the tight-binding representation with attached source or detector 
terminal is introduced. The formalisms for calculating the local particle density, the local 
current density vector and the local conductance are then derived. Sec. Ill gives numerical 
results along with some discussions. Finally, a brief summary is presented in Sec. IV. 
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II. MODEL AND FORMALISM 



In order to study the scattering due to PNJ, we consider two kinds of open bipolar 
graphene systems (armchair and zigzag ribbons) as shown in Fig.l. The bipolar graphene 
ribbon consists of semi-infinite electronlike ribbon [orange region] and semi-infinite holelike 
ribbon [green region] along x-direction with a sharp p-n interface located at x = 0. Electron 
flow is injected into graphene system from a source lead located at (— L, 0) in the n-region. 
Here we assume that the source lead and the bipolar ribbon are in contact with six lattices 
[see the blue area in Fig.l]. The injected electrons in the n-region can spread in all directions. 
Because of the open boundary condition, left-going electrons can finally escape into infinite 
graphene electrode while right-going electrons [shown in Fig. 2] can then be scattered only 
by p-n interface [thick black line]. Consequently, the response signals are converged around 
the symmetric site (L, 0) [red area in Fig.l]. In order to investigate the focusing current, we 
couple a detecting electrode locally in the p-region and study the local current (conductance) 
flowing from that electrode. Clearly, the local current depends on the coupling position of 
the detecting electrode. 

The total Hamiltonian including the infinite graphene ribbon in the tight-binding 
representation^ and the source/drain electrode that is expressed in k space with the free 
electron model can be written as: 



where i = (i x , i y ) is the index of the discrete site on the honeycomb lattice which is sketched 
in the Fig.l, and a\ and a\ are the annihilation and creation operators at the site i. Here e\ in 
the first term of Eq.(l) is the on-site energy (i.e., the energy of the Dirac point) which can be 
controlled experimentally by the gate voltage. In the n-region or p-region far away from PNJ 
all the on-site energy are the same with ei = E n or ei = E p . Near PNJ, e\ changes from E n to 
E p abruptly or smoothly for the sharp PNJ or smooth PNJ, respectively. The second term in 
Eq.(l) is the nearest neighbor hopping term with the hopping energy t, '< ij >' denotes the 
nearest neighbor lattice sites. When the graphene ribbon is under a uniform perpendicular 
magnetic field B z = B, a phase </>y is added in the hopping term, and </>y = J1 J A ■ dl/(fio 
with the vector potential A = (—By, 0,0) and the flux quanta O = fi/e. Finally, the last 




<ij> 




(1) 
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term in Eq.(l) represents the Hamiltonian of the source and detector leads described in the 
k space and their coupling to the graphene lattices i a . Here a = s,d represent source and 
detecting electrodes and d a ^ (d a k ) is the annihilation (creation) operator of the electrons 
in the electrode a. 

When the electron flow is injected from the source electrode into the graphene in the n- 
region, the response signal is induced everywhere in the p-region. To make a thorough study 
on the focusing effect, we consider three physical quantities in the p-region: (1) the local 
current density vector, (2) the local particle density, and (3) the local current (conductance) 
of the detecting electrode. For quantities (1) and (2), we consider the system without the 
coupling of the detecting electrode so that the influence of the detecting electrode can be 
eliminated. 



A. Local current density vector 

The general current density vector Jy from the site i to its nearest neighbor site j can be 
expressed as:— 

h = \\ dE [G^^Hji-HyGfj^)] 

= 2^Re J dE [te^Gf^E)] (2) 

where e is the electron charge, G^ is the matrix element of the lesser Green's function of the 
scattering region. Because the graphene ribbon is translation invariant in the p/n-region, 
the central scattering region can be chosen arbitrarily as long as the source sites \ s and the 
detector sites i<j are included. From the Keldysh equation, the lesser Green's function is 
related to the retarded and advanced Green's functions, 



G < (E) = G r (E) 



G a (E) (3) 



Here the sum index a = L,R,s denote the left/right graphene lead and source lead with 
a 7^ d because of the decoupling of the detector lead. The retarded Green's function G'\E) = 
[G a (-E')]^ = {El — H — Yla ^aC^)} _1 > wnere H is Hamiltonian matrix of the central 
scattering region and I is the unit matrix with the same dimension as that of Ho, S„ is 
the retarded self-energy function from the lead-a. 17 a can be obtained from 'E r L , R (E) = 
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H-c,L/Rg r L/R (E)H L/RjC , where H c>L/jR (H L/H)C ) is the coupling from central region (lead-L/R) 
to lead-L/R (central region) and g r L ^ R (E) is the surface retarded Green's function of the 
semi-infinite lead which can be calculated using a transfer matrix method.— Concerning the 
source lead, we take the wide band approximation, then the non-zero elements of self energy 
matrix Ti r s (E) = —iT s /2 is energy independent, where line-width function T s = lurf p s (Ep). 

in Eq.(3) is the lesser self energy of the lead-a. Because the isolated lead is in the 
equilibrium, can be obtained from the fluctuation-dissipation theorem: 



K(E) = [K(E) ~ K(E)} fa(E) 
= iT a (E)f a {E) 



(4) 



with H a a = Xftt and f a (E) = f (E - eV a ) where f (E) = l/[exp(£ /k B T) + 1] is the Fermi 
distribution function. V a is the external bias in the terminal-a. Since we are interested 
only in the local response due to the source flow, the external biases are set as V s = 5V 
and Vl/r = 0. In calculating transport properties, we divide G < (E) into equilibrium and 
non-equilibrium parts as 



G < {E) = G r (E) 



+ G r (E) 



if (E)J2^(E) 



G a (E) 



j2{U E )-fo(mr a (E) 



G a (E) 



(5) 



where the equilibrium term does not contribute to the transport and can be dropped out 
from now on. It is the non-equilibrium term that gives rise to the system response to the 
electron injection from the source lead. Because of V s = 5V and Vl/r = 0, we have 



G<(E) = iG r (E) [f s (E) - f (E)} T s G a (E) 



(6) 



Substituting Eq.(6) into Eq.(2) and considering the limit of small source bias, the local 
current density vector Jy (or the local conductance density vector 3ij/5V) can be expressed 
in the following form at zero temperature: 

2e 



h/5V = -^Imjte^ [G T {E F )iT s G a {E F )\.^ 



(7) 



It should be noted that the current density Jy in Eq.((7j) is defined between the lattice sites 
i and j with the direction from site i to j. In order to obtain the local current density vector 
J ; at the site i, we take the weighted average on Jy over all the nearest neighbors j. 



B. Local particle density 



The local particle density (i.e. the electron occupation number) is defined as 



Pi = ~ 




(8) 



where Gj? is the diagonal element of the lesser Green's function G < in Eq.(3). Similar to 
the derivation of local current density vector J^/SV, here we consider only the variation 
of the local particle density caused by the electron injection from the source lead. At zero 
temperature and the small bias 5V limit, substituting Eq.(6) into Eq.(8), the variation of 
the local particle density is expressed as: 



Since the Hamiltonian is defined at discrete lattice sites, the local quantities can also be 
defined at each lattice site. Such a local quantity is feasible but not necessary. In fact, for 
graphene, we can define the "local" quantity by averaging over six discrete sites in a unit 
cell of honeycomb lattice. This average can eliminate the strong variation of local quantities 
in the A and B sublattices. Now every local site can be determined from the coordinates 
(x, y) shown in Fig.l. For example, the local injection area displayed in a blue area in Fig.l 
is located at (—3, 0) in Fig. 1(a) and (—5, 0) in Fig. 1(b). In the whole lattice region in Fig.l, 
there are 7x3 and 11x3 units in Fig. 1(a) and Fig. 1(b), respectively. 

C. Local conductance 

Concerning the local conductance, the detecting lead-d is coupled to the graphene ribbon 
in the p-region. Similar to the local particle density, here the detecting lead also couples to 
six sites i d in a unit cell of honeycomb lattice. The current flowing to the detecting lead-d 
can be expressed as 



S Pi /5V = [ Pi (V s = SV) - Pi (V s = 0)]/5V 
= i- [G r (E F )T s G a (E F )] ii 



(9) 



tl=*2 




(10) 
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Using the Dyson equation in the time contour, we can get the Landauer-Biittiker formula^ 
which is expressed in terms of non-equilibrium Green's functions: 

Jd = TXj ^ T *AEM(E) - fa(E)]. (11) 

a 

with a = L, R, s representing the left/right graphene leads and source leads. Since we shall 
concentrate only on the response current induced by the current injected from the source 
lead, we use the following boundary conditions Vl,r = Vd = and V s = 5V. The current 
now becomes, 

J *=lJ ^[T d AE)(f d (E) - fs(E))\. (12) 

where T^ s is the transmission coefficient from the source lead located at the site i s 
to the detecting lead located at the site which can be calculated from Td tS (E) = 
Tr [r,iG r (E)r s G a (E)], where G a = G r * is the advanced Green's function in the scatter- 
ing region. In the wide band limit, the line-width function T s / d (E) = i{^ r s / d — ^TL) = 
27T7 2 'p s /d(E F ). Here f s /d(E) in Eq.(12) is the Fermi distribution function of the source and 
detecting lead, and f s (E) = f (E — eV s ) and fd(E) = fo(E). Considering the zero tempera- 
ture and small bias V s limits, local conductance contributed by the source electron flow can 
be expressed as: 

G id = J d /5V = jT d , s (E F ). (13) 

III. NUMERICAL RESULTS AND DISCUSSION 

In the numerical calculations, we set the nearest-neighbor carbon-carbon distance a = 
0.142nm and the second nearest-neighbor distance b = \/3a ~ 0.25nm, the hopping energy 
t = 2.75eV as in a real graphene sample.- In this paper, we consider only the focusing effect 
of the symmetric PNJ, in which the electron density of n region is the same as the hole 
density in p region, i.e., p e = ph- For simplicity, we set Ep = 0. Hence in the n-region far 
away from PNJ the on site energy ti = E n = — E while e- x = E p = E in the p-region far 
away from PNJ. Near PNJ, e- Y changes from — E to E abruptly (smoothly) for the sharp 
(smooth) PNJ. 

Experimentally, it is more convenient to measure the electric conductance. So, in order 
to detect the focusing effect by a single PNJ in graphene, one can use a small electric 
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contact (such as a STM probe) as a source of electron flow in the n-region and another 
local probe located in the p-region as a detector. Electric conductance between the two 
contacts measures the transmission probability for a charged carrier from the source to the 
detector. Numerically, We have calculated the local conductance G- ld and confirmed that 
the distribution of local conductance is similar to that of the local particle density 8p\j8V . 
For this reason, only the numerical results on local particle density are shown in this paper. 
In addition, in order to visualize the focusing process, we also show the distribution of local 
current density vector in the p-region. 



A. Focusing effect in very low energy regime 

Now we study the focusing effect in the graphene ribbon with a sharp PNJ. For a zigzag 
ribbon or a armchair ribbon with sharp and symmetric PNJ, the spacial distribution of the 
local particle density Sp(x,y)/5V in p-region due to electrons coming from the source lead 
is shown in Fig.3(a) and Fig. 3(b), respectively. Following observations are in order. First 
of all, electrons injected at (— L, 0) in n-region can be focused around (L, 0) shown as red 



spot in Fig. 3 which is similar to Ref. 13(. This is because the Fermi energy Ep is close to 
Dirac energy E = 0.05t so that the energy dispersion is nearly linear, i.e., E ~ kb^-t 
where k is module of momentum vector k. The charged carriers scattering through PNJ 
can mimic the refraction of light by left-handed metamaterials with refraction index equal 
to -1. Secondly, besides the focusing spot (red and green region), there is also a weak 



interference pattern (blue wave pattern) shown in Fig. 3, which is different from Ref. 13( in 
which the wave pattern is absent. In fact, the wave pattern is solely due to the boundary 
of the nanoribbon. When an electron is injected from the source area, it can propagate 
in all directions and the right-going electrons can be scattered by either the boundary of 
nanoribbon (thin black lines) or sharp PNJ (thick black line) as shown in Fig.2. In the 
p-region, interference pattern is due to the interference between the state scattered by both 
boundary and PNJ (thick blue lines) and the state scattered only by PNJ (thick red lines). 
The spacial period of the interference is proportional to the momentum k or E . Finally, the 
focusing phenomena in zigzag ribbon is slightly different from that of armchair ribbon: the 
electron flow is perfectly focused in armchair ribbon [panel (b)], but can't be fully focused 
in the zigzag ribbon [panel (a)]. This is because the energy band structures are different 
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for the armchair ribbon and zigzag ribbon. In the following, we will examine the different 
focusing effects in detail. 

B. Focusing in zigzag ribbon 

When Fermi energy is gradually moved away from Dirac point, the energy spectrum is 
not ideal conical anymore. In Fig.4 we plot the contour lines of dispersion relation E(k x , k y ) 
of graphene sheet^Z with energy interval between nearest contour lines 5E = O.lt. The 
panel (a) is for the graphene sheet with the carbon-carbon bond along the x-direction which 
corresponds to the armchair graphene ribbon, and the panel (b) is for the graphene sheet 
with the carbon-carbon bond along the ^-direction corresponding to the zigzag ribbon. 
The deviation of ideal conical energy spectrum is clearly exhibited even at small energy 
E = 0.2t (the second small contour lines around the Dirac points K and K' show anisotropy 
behaviors). Since the p-n interface at x = is along y direction, the y component of 
momentum, k y , is conserved during the scattering. As a result, the incident wave vector 
k Xt i n , the reflecting wave vector k Xyr , and the transmitting wave vector k x j must lie on the 
black dotted lines in Fig.4. When an electron with energy E = Ep, velocity (v x ,v y ) and 
corresponding momentum (k x ^ n , k y ^ n ) with respect to Dirac point K, injects from n-region 
and is scattered at the p-n interface, according to the identical direction of V x , we can solve 
the reflecting and scattering momentum k x>r and k X)t using the energy conservation and k y 
conservation. For the zigzag ribbon [corresponding to Fig.4(b)], k Xtin can be intra-scattered 
to k x ^./ t around K [K = (1^,0)] valley or inter-scattered to k x ^ r / t in K' [K' = — (f^-,0)] 
valley. The interband scattering states is symmetric which satisfies k' x r , t = —k x ^ n . The 
intraband scattering, however, exhibits asymmetric properties that 8k x ^ r / t = k Xj i n +k Xjr u ^ 
or v x> i n 7^ — v Xir , v Xt i n ^ for any fixed k y . In Fig.5, we plot the asymmetric relation 
between k x ^ n and k l ^ r ^ in the intraband scattering and the symmetric relation between 
k x ^ n and &* n *% in the interband scattering. We see that in the intraband scattering, the 
larger k y , the larger derivation is, while in the interband scattering, k l ^J r t is always 

equal to —k Xtin for all k y . It is known that interband scattering is weak— in pure samples 
due to the large momentum shift, so the asymmetric intraband scattering is dominant in 
zigzag ribbon PNJ. As a result, the refraction index can not be strictly equal to —1 and 
the charge flow can't be fully converged at the symmetric spot. In Figj6] focusing effect for 
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E = O.lt and Eq = 0.2t are plotted, respectively. Comparing Figj6]^a) and Figj6]^b) we see 
that it is more difficult to focus the electron beam for larger momentum k (energy Eq). 

C. Focusing in armchair ribbon 

For the armchair ribbon [corresponding to Fig.4(a)] only intraband scattering occurs. 
Now k Xj i n can be symmetrically intra-scattered to k x r /t in K valley with —k Xt i n = k x>r = 
k X)t since k x in Fig.4(a) is symmetric about k x = 0. So, although the energy dispersion 
of armchair ribbon is also not strictly linear at high energies as in zigzag ribbon, due to 
the symmetric scattering, the focusing effect is always perfect in armchair ribbon for all 
Eq < t (i.e. in Dirac fermion regime). Furthermore, with increasing of Eq, the electron 
flow coming from n-region shows better convergence in the p-region with smaller focusing 
spot and stronger intensity. Furthermore, The spacial period of the interference pattern is 
proportional to the momentum k or Eq, which can be clearly seen by comparing Fig. 3(b), 
FigfT^a) and Figfftb). . 

Roughly speaking, two energy regimes are considered for the Dirac Fermion according 
to band structure of graphene: (1). 'Near linear dispersion' regime < Eq < 0.5£ where 
Eq rs kb^t. (2) 'Beyond linear dispersion' regime 0.5t < E < t where the energy spectrum 
is non-conical. The focusing effect corresponding to these two regimes are plotted in FigfT] 
and FigJHl respectively. In the first regime, with the near linear dispersion relation, velocity 
v x or Vy is roughly a constant, and ky^ n j k x ^ n ~ Vy^ n j v X i n , kyi /k x>t « Vy tt /v Xit . For the 
symmetric scattering (k x> i n = k X) t) in the armchair ribbon, refraction index n ~ — 1 giving 
rise to the convergent spot shown in FigJTl When Eq is large enough (larger than 0.5t), 
energy spectrum is non-conical and velocity now depends on momentum. This leads to a 
different focusing effect shown in Figj8]in which crossed focusing zone is present. 

In order to show the focusing of the electron flow vividly, instead of the contour of 
local particle density in FigJT^b) and FigJH^b), the quiver of local current density vector 
around the convergence spot in p-region are plotted in Figj9] and FigflO| respectively. For 
demonstration purpose, the local current density is plotted at every other site. The arrow on 
each site denotes the local current density vector whose module and direction are described 
by the size/color and orientation of the arrow respectively. In the low energy regime (FigJHl), 
the vectors of local current density converge conically to the focusing spot [red spot in 
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FigJT^b)]. On the other hand, current density is converged mainly from four crossed corners 
in the high energy regime (Fig JTOj) . Furthermore, comparing Figj9] and FigJTUl it is clear 
that electron flow with larger E gives better convergence. 

D. Effect of disorders in armchair nanoribbon 

As discussed in the previous sections, the clean graphene PNJ is investigated. In a real 
device, the disorder is always present. In this subsection, we study the disorder effect on 
the perfect focusing in the armchair nanoribbon. We consider two kinds of disorders, one is 
induced by random on-site potential 5ei and the other is due to the edge defect.-^ The random 
on-site potentials 8e\ with a uniform distribution [— w/2, w/2] are added near PNJ within 
the width of 18a where w is disorder strength. The edge defect is modeled through missing 
atoms on the graphene edge. We model the missing atom by setting the corresponding 
hopping matrix elements to zero. The edge roughness is controlled by p, the probability of 
a missing atom on the outermost row [the red line in Fig. 11]. For both on-site potential 
disorder and edge defect, all data are obtained by averaging over 500 configurations. 

In Fig{12] we plot the contour of local particle density in armchair ribbon with a sharp 
PNJ for Eq = 0.5t [same as in FigJH](a)] in the presence of random on-site potential disorder. 
In FigJT2"| the width of armchair ribbon is set to 1056, the source flow is injected from the 
honeycomb unit cell at (— 210a, 0) and focused around the spot located at (210a, 0). Panel 
(a), (b), (c) and (d) correspond to different disorder strengths w = 0, 0.2, 0.5 and 1.0. For the 
small random potential strength w (e.g. w = 0.2), the interference pattern and the focusing 
spot can be well kept. On the other hand, for the large w (e.g. w = 1.0), we can see that 
the random potential disturbs the interference between forward and backward scattering, so 
the interference pattern is smeared, which increases the density of state outside the focusing 
spot. Consequently, the intensity of focusing spot decreases. However, we emphasize that 
although random potential disturbs the interference pattern and reduces the intensity of 
the focusing spot, the focused spot is clearly visible and its size still remains unchanged. It 
means that the focusing effect is robust against random potential, especially in the weak 
disorder case. 

In FigHni we plot the contour of local particle density in the presence of edge defect 
for Eq = 0.5t. Panel (a), (b), (c) and (d) are corresponding to the different probability 
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p of a missing atom on the outermost row with p = 0, 0.1, 0.2 and 0.5. We find that in 
the presence of edge disorder, the size of focusing spot increases clearly and the focusing 
intensity is greatly reduced. So the effect of the edge defect on the focusing effect is more 
significant than that of the random potential. But the focusing spot and interference pattern 
still survive and are clearly visible (see FigJTBTi). even in the strong edge defect case with 
p = 0.5. 

To estimate the disorder strength needed to reduce the intensity of focusing spot, in 
Fig lMT a) and (b), we plot the maximum value (the value at the focusing spot central [L, 0]) 
of the focused spot LDOS, mai . vs strength of random potential w and the probability of a 
missing atom p. Considering the computational cost, here we take 200 configurations and 
label the error bar. From FigJH^a), we find that for weak random potential (when w < 0.5), 
LDOS max hardly changes with w, and focusing effect remains unchanged [see Fig JTW a). (b) 
and (c)]. Beyond the weak disorder regime (w > 0.5), LDOS max declines abruptly, and 
focusing effect can't be kept as good as in the weak disorder regime [see Fig fl27 d)]. On 
the other hand for the edge defect (see Fig jMf b)). we can see that the electron beam can 
be focused perfectly at p = and p = 1 because the graphene ribbon edges are intact at 
both p = and 1. When p increases from to 1, more and more atoms in edge are missing 
until two edges are completely peeled. Correspondingly LDOS max decreases first and then 
increases since the edges are the most random when p is around 0.5. We notice in Fig fUI b). 
comparing to LDOS max near p — 1, LDOS max is reduced faster near p = 0. It means that the 
vacancy defect (a few atoms are missing on edges) destroys focusing effect more significantly 
than the adsorption defect (a few atoms are attached to edges). 

E. Focusing of armchair ribbon with smooth PNJ 

Up to now, we have studied focusing effect by the sharp PNJ. But in realistic graphene 
based PNJ or PNPJ, the potential changes smoothly from E n to E p within a width d. 
The width d is of the order of the separation between the graphene layer and the top gate 
and d ~ tens nm.- In such a smooth PNJ, backscattering is present near PNJ in the 
distance proportional to k y and interface width d, which reduces the possibility of Klein 
tunneling. For a linear electrostatic potential U(x) = (vkp/d)x, the angular dependent 
transport probability 11 T{6) = e -^( k Fd) sm 2 (e) ^gj-g q j s incident angle. It is obvious that 
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the smooth PNJ will reduce the intensity of the focused electron beam due to the decreased 
transport probability T(9). It appears that it also increases the size of the focused spot. 
Actually, it is not the case due to the following reason. For a single n-p junction (whether 
smooth or sharp), the electrons (holes) with an energy equal to the chemical potential \x = 
and momentum k x = hp cos(6 l ) transport from conduction (valence) band to the the valence 
(conduction) band with the conserved k y = kpcos(9) but k' x = —k x , leading to the almost 
unchanged focusing pattern, as shown in Figfl5l 

The smooth PNJ can be modeled by smoothly varied Dirac energy U(x) across the PNJ. 
In the numerical calculation, we use the following U(x), 



U(x) 



-E [1 + sinh(f )/sinh(£=fa)] , x < 

(14) 

E [1 - sinh(^)/sinh(2±p)] , x > 



where L is the distance between source probe located at (— L, 0) and PNJ at x = 0. In 
Fig JTBT a). with L = 320 x 0.43nm, U(x) for different xo have been plotted. In Fig JToT b. c, 
d) with the same parameters used in FigJ7](a) in which the sharp PNJ is used, the contour of 
local particle density is re-plotted for the smooth PNJ shown in Fig JToT a). We can see that 
the quantum interference between the states scattered by forward and backward scattering 
is present. The density of state outside focused spot is increased. Moreover, intensity of the 
convergent spot is reduced comparing to that of sharp PNJ due to the reduced transmission 
probability. The wider the PNJ interface, the more the focusing effect is reduced because of 
the smaller T(8). For example, when the PNJ width d = 0, the maximum of local particle 
density of state LDOS mai = 0.0095 [see FigJ7(a)], increasing d gradually, LDOS max = 0.0068, 
0.0035 [see Fig JTBT b) and (c)]. When the PNJ width d is reaching to Fermi wavelength [such 
that k F d ~ 1, in Figjl5l kp ~ l/(15a)], the intensity of focused electron beam decrease very 
slowly [see Fig JToT c) and (d)]. For the very big d » 1/kp, the intensity of focused spot 
decreases and its size increases continually. In addition, due to the Klein tunneling, the 
focusing effect can still occur and the convergent contour is almost the same as that of the 
sharp PNJ, although the intensity of focusing spot is reduced. 
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F. Focusing of armchair ribbon in the presence of small perpendicular magnetic 
field 

In the presence of small perpendicular magnetic field, the momentum k y is not a con- 
served quantity. In this case, electrons and holes are deflected in opposite directions due 
to the opposite Lorentz force. So the injecting electron flow in the n-region now can't be 
effectively converged in the p-region and the focusing spot is smeared. Considering the op- 
posite deflection for the electrons and holes, the smeared convergent spot will move away 
from symmetric point (L, 0). The larger the size of scattering region, the focusing effect 
is less significant because of the more deflection in the larger size. On the other hand, 
the transmission probability T of a single PNJ becomes magnetic field dependent on the 
field scale B* = (h/e)kp/d with which the cyclotron radius l cyc i = hkp/eB becomes com- 
parable to the width d of PNJ. The maximum angle rotating away from normal incidence 
9-max — i arcsin (5 /£?*). The transmission probability of a bipolar ribbon is suppressed as 
T(B <B*)(x^(l- (B/B*) 2 ) 3 / 4 where d is width of PNJ and W is width of ribbon. 31 The 
influence of transmission probability on focusing effect is mainly to reduce the intensity of 
focused spot. So in the presence of the weak magnetic field, not only the intensity of focused 
spot is reduced but the focusing pattern is destroyed by the deflection of electron beam as 
well. 

In FigJT6| we plot the local particle density in the armchair ribbon with small magnetic 
field for sharp and symmetric PNJ. For the sharp PNJ, intensity of focused spot is reduced 
not so severely as in the smooth PNJ. The magnetic field B is expressed in terms of magnetic 
flux BSq in the unit of (fio/n where So = fv^ 2 is the area of a honeycomb unit cell and 
0o = V e i s the flux quanta. Here BSq = O.OOOl^o/ 71 " corresponds to the magnetic field 
B = OAT. From FigJT6| we can see that in the low energy regime [panel (a), Eq = 0.2t] 
the focusing effect is destroyed severely, and focusing effect is hardly destroyed in the higher 
energy regime [panel (b), E = 0.9t] due to the less influence of magnetic field on electron 
flow with higher energy. 
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IV. CONCLUSION 



In conclusion, using the tight-binding Hamiltonian, we report the focusing of electron 
flow in zigzag or armchair graphene ribbon with a single symmetric PNJ. For a sharp PNJ, 
in the very low energy regime (\E F — E n / p \ = E < 0.05), graphene ribbon exhibits almost 
conical energy spectrum and the electron flow coming from n-region can be converged in 
the p-region perfectly. When energy E increases, however, energy spectrum gradually 
deviates from linear behavior. And the band structures are different for zigzag ribbon and 
armchair ribbon. For the zigzag ribbon, although the interband scattering is symmetric 
with —k l x n ^ = k™lj t but the dominant scattering, intraband scattering, is asymmetric with 
k™in a + ^Tr V /t 7^ 0- As a result, the electron flow coming from source in n-region can't 
be converged in the p-region. As for the armchair ribbon, only the intraband scattering 
exists which is always symmetric with —k x ^ n = k x ^ r / t for all k y . This leads to a perfect 
focusing effect for all energy Eq < t regardless of linear or nonlinear dispersion relationship 
of Dirac Fermion. Specifically, the electron flow converges conically in the low energy regime 
(E < 0.5t), and converge mainly from four crossed corners in the high energy regime 
(0.5t < Eq < t). When disorder is present, the perfect focusing in the bipolar ribbon can be 
robust against disorder induced by the random potential. The perfect focusing is however 
drastically affected by the edge defect, not only the intensity of focused spot is reduced, 
the size of spot is also increased. Furthermore, when the real smooth PNJ is considered, 
Klein tunneling is reduced significantly due to the backscattering. In this case the intensity 
of convergent spot is reduced, but the convergent contour still remains the same. Finally, 
small perpendicular magnetic field deflects the electrons and holes in opposite directions, 
which destroys the perfect focusing effect especially in the low energy regime. 
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FIG. 1: (Color online) A schematic diagram of the graphene PNJ in a armchair ribbon [panel 
(a)] and in a zigzag ribbon [panel (b)]. The graphene ribbon is along x-direction and sharp p-n 
interface is located at x = 0. The electron flow injected at (— L, 0) [the blue area] in the n-region 
[orange lattice region] is focused around the symmetric site (L, 0) (the red area) in the p-region 
[green lattice region]. 
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FIG. 2: (Color online) After injecting from the source lead in the n-region, the right-going electrons 
are scattered by the PNJ and boundary of the bipolar graphene ribbon. 
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x (b=0,25nm) x (3a=0.43nm) 

FIG. 3: (Color online) Distribution of local particle density Sp(x, y)/5V in a graphene ribbon with 
a single sharp PNJ at x = 0. Panel (a): zigzag ribbon with ribbon width W = 500 x 3a. The 
source flow is injected from the honeycomb unit cell at (—8006,0) and focused around the spot 
located at (8006, 0). Panel (b): the armchair ribbon with ribbon width W = 5016. The source flow 
is injected from the honeycomb unit cell at (—320 x 3a, 0). The other parameter used: Eq = 0.05i. 
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FIG. 4: (Color online) The contour of dispersion relation E(k x , k y ) of graphene sheet. The energy 
interval between nearest contour lines is 5E = O.li. Panel (a): E(k x , k y ) of the graphene sheet with 
the carbon-carbon bond is along the x-direction corresponding to armchair ribbon or Fig.?? (a). 
Panel (b): E(k x ,k y ) of the graphene sheet with the carbon-carbon bond is along the y-direction 
corresponding to zigzag ribbon or Fig. 1(b). 
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FIG. 5: (Color online) Scattering momentum k 1 ™^,™ and &* n *y[ vs injecting momentum k x ^ n for 



the zigzag ribbon with sharp PNJ. In the interband case, 



-k x ,in for all conserved k y (red 



line). While in the intraband case (black lines), k l ^y t is not equal to —k x ^ n . Different black lines 
along the black arrow are corresponding to k y = 0, 0.1/6, 0.15/6, 0.2/6, respectively. 
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x (3a=0.43nm) x (3a=0.43nm) 

FIG. 7: (Color online) Contour of local particle density in armchair ribbon with a sharp PNJ for 
E Q = O.lt [panel (a)] and E = 0.2t [panel (b)]. 
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x (3a=0.42nm) x (3a=0.43nm) 

FIG. 8: (Color online) Contour of local particle density in armchair ribbon with a sharp PNJ for 
E Q = 0.5* [panel (a)] and E = 0.9* [panel (b)]. 
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FIG. 9: (Color online) Instead of contour of local particle density in Figj7]Jb), the quiver of local 
current density vector around convergence spot is plotted. 
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FIG. 10: (Color online) Instead of contour of local particle density in FigJHJ^b), the quiver of local 
current density vector around convergence spot is plotted. 
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FIG. 12: (Color Online) Contour of local particle density in the armchair ribbon with a sharp PNJ 
for Eq = O.bt. Panel (a), panel (b), panel (c) and panel (d) are corresponding to random on-site 
potential strengths w = 0, 0.2, 0.5 and 1.0, respectively. 
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FIG. 13: (Color Online) Contour of local particle density in armchair ribbon with a sharp PNJ for 
Eq = 0.5t. Panel (a), panel (b), panel (c) and panel (d) are corresponding to p = 0, 0.1, 0.2 and 
0.5, respectively. 
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FIG. 14: According to Figll2l and Figll3|, maximum value (the value at the focusing spot central 
[L, 0]) of the focusing spot LDOS maa; vs the strength of random potential w [in panel (a)] and 
probability of one missing atom p [in panel (b)] are plotted, respectively. 
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FIG. 15: (Color Online) panel (a): smoothly changed U(x) forming smooth PNJ used in panel (b, 
c, d), in which xo = 0.25 * 3a, 0.75 * 3a and 3 * 3a, respectively. Panel (b, c, d): Contour of local 
particle density in armchair ribbon for Eq = 0.1 [same to FigJT^a) where sharp PNJ is used]. For 
different panels, smooth PNJ shown in panel (a) are used respectively. 
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FIG. 16: (Color online) Contour of local particle density in the p-region for a armchair ribbon in 
the perpendicular weak magnetic filed. The sharp PNJ is located at x = 0. The magnetic field 
BS = O.OOOl^o/vr, E = 0.2t in panel (a) and E = 0.9* in panel (b). 
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